Untwisting of a cholesteric elastomer by a mechanical field 
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A mechanical strain field applied to a monodomain cholesteric elastomer will unwind the helical 
director distribution. There is an analogy with the classical problem of an electric field applied to 
a cholesteric liquid crystal, but with important differences. Frank elasticity is of minor importance 
unless the gel is very weak. The interplay is between director anchoring to the rubber elastic matrix 
and the external mechanical field. Stretching perpendicular to the helix axis induces the uniform 
unwound state via the elimination of sharp, pinned twist walls above a critical strain. Unwinding 
through conical director states occurs when the elastomer is stretched along the helical axis. 



PACS numbers: 61.30.-v, 61.41. 
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Monodomain cholesteric elastomers are formed by 
crosslinking mesogenic chiral polymers in the cholesteric 
state with a properly formed helical director twist. The 
subsequent retention of the helical state as an elastic 
equilibrium Q is a natural consequence of topological 
imprinting of textures in the crosslinked network, seen 
in a number of other elastomers with liquid crystalline 
order and other microstructure. Recently an interesting 
aspect of chiral imprinting was established by crosslink- 
ing nematic polymers in a chiral state purely induced by 
a chiral solvent B]; on removal of the solvent, the network 
of chemically achiral nematic chains remains macroscop- 
ically cholesteric. Such an imprinting has been envisaged 
long time ago |s| on phenomenological grounds. It is now 
important to consider the mechanical possibilities of such 
solids with a helical microstructure, expecting new tran- 
sitions and instabilities characteristic of liquid crystalline 
elastomers. Additionally there are obvious device appli- 
cations of such materials, which combine all the optical 
properties of twisted nematic liquids with the remarkable 
mechanical characteristics of rubbers. There is some ex- 
perimental evidence that such effects are indeed ob- 
servable and our hope is that this theoretical work will 
stimulate more studies in this field. 
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FIG. 1. The initial director 110(2) in a cholesteric helix 
makes an azimuthal angle <f>o = qoz with the x-axis; the heli- 
cal pitch is p — ir/qo. Two principal directions of mechanical 
deformation, A xx and A zz , are shown by arrows. 



Consider a monodomain cholesteric elastomer with an 
ideal helically twisted director no (z) in the xy plane, ini- 



tially making angle <fio — qqz with the x axis, Fig. @. We 
shall examine two specific cases of imposed uniaxial ex- 
tension: (i) the transverse deformation A xx = A, in the 
plane including no, and (ii) the longitudinal deformation 
along the helix axis A zz = A. 

The symmetry obvious from Fig. [j] requires that in 
the case (i) the director remains in the xy plane, char- 
acterised by the azimuthal angle <f>(z), while in the case 
(ii) one may expect a conical texture with n(z) inclined 
towards the stretching axis z and, therefore, described 
by two angles 8 and <f> (cf. Fig. || below). In ordinary 
liquid cholesterics subjected to, e.g., a magnetic field H z , 
such conical states are not generally seen, preempted by 
the 90°-switching of the helix axis and then untwisting in 
the "transverse" geometry ||. We shall see that in elas- 
tomers, due to the chiral imprinting, this regime is not 
possible and the conical director configurations should 
occur. 

An elastic material with a microstructure represented 
by an independently mobile director orientation is anal- 
ogous to the Cosserat medium. In the limit of linear 
elasticity the relative rotation coupling between the di- 
rector rotation ui = [n x 5n] and the antisymmetric part 
of strain, fli — eyk^jk, 
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D x [n x (SI - u>)] 2 + D 2 n ■ e (s) ■ [n x (SI - w)], (1) 



has been first written down phenomenologically by de 
Gennes being the symmetric part of the small 

strain defined as e = A — 6. This symmetry-based ex- 
pression is only valid for small deformations, having only 
linear and quadratic terms in the local relative rotation. 

The microscopic statistical-mechanical theory of ne- 
matic rubber elasticity, e.g. 0, obtains a generalisation 
of the classical rubber-elastic energy density in the form 
of a complete frame-independent expression 



F =^Tr(| o -A T -r 1 -A 



(2) 



plus the constraint of material incompressibility, ex- 
pressed by the condition Det (A) = 1 on the strain 
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tensor. Apart from the strain tensor, the other en- 
tries in the Eq. (§) are £ Q = £±S + (£|| — ^jJn no and 

r 1 = (l/^x)|+(l/^||-l/^x)nn, the reduced shape and 
inverse shape tensors characterising the Gaussian distri- 
bution of nematic polymer chains before and after the 
distortion A. The rubber shear modulus fi = n s kBT 
(with n s the number density of network strands, propor- 
tional to the crosslink density) is that characteristic of 
the underlying isotropic rubber and sets the energy scale 
of distortions. The free energy density Eq. (||) is known 
to be valid up to large strains and correctly predicts the 
opto-mechanical responses and the soft elasticity of ne- 
matic elastomers. The free energy F is a function only of 
the chain anisotropy r = £\\/£±, the ratio of the effective 
step lengths parallel and perpendicular to the director, ft 
is an independently measured parameter accessible from 
neutron scattering or from spontaneous mechanical dis- 
tortions on going from the nematic to isotropic phase. 
Unless there is a large nematic order change induced by 
A, the shape £ is essentially just a rotated version of £ Q , 
a uniaxial ellipsoid with the long axis (at r > 1) oriented 
along n instead of no. 

Embedded in the general expression Eq. (g) is the 
penalty for local director deviations from the orientation 
no imprinted into the network at formation. When no 
elastic strains are allowed, this elastic energy reduces to 

3 1 (r- l) 2 , 
F - ^ + sin 9 ( 3 ) 

where is the local angle between n and no. The elastic 
penalty for such a deviation, appropriately proportional 
to the square of chain anisotropy, is the coefficient Di 
of the de Gennes' phenomenological expression at small 
deformations, Eq. ([!]). This has to be compared with the 
Frank elastic penalty for director curvature deformations, 
|X(Vn) 2 . The length scale £ ~ y/K/fi at which the 
two energy contributions are comparable is usually small: 
£ - I(T 8 m for a typical K ~ fO- n J/m, fi ~ f0 5 J/m 3 
and not too small anisotropy, r. This is rather less than 
the cholesteric pitch p, which is a characteristic scale in 
our problem. Therefore, the anchoring of the director 
n to the rubbery matrix, described by Eq. (0), tends to 
dominate over Frank effects. 

We shall assume that a cholesteric elastomer is locally 
like a nematic in its elastic response: rubber elasticity 
is determined on the scale of network crosslink separa- 
tions (a few nanometers) , whereas cholesteric pitches are 
10 3 times longer. We can at once see why the chiral 
structure is stable and how mechanical fields can desta- 
bilise it. With no elastic strain, the free energy penalty 
is ~ hD±(<p — 4>o) 2 for rotating the director away from 
its original helical texture 4>o — Qqz. On the other hand, 
if strains are applied, the rubber can lower its elastic en- 
ergy Eq. (|^) by rotating the director n towards the axis 
of principal extension. This general principle of adjust- 
ing the microstructure to minimise the clastic energy is 
seen in its ultimate form in the effect of soft elasticity 



when a stretched nematic rubber may reduce its 
effective modulus (the slope of a stress-strain curve) to 
zero by optimising the director rotation and associated 
shear strains. 

Distortions in a cholesteric elastomer cannot be soft, 
because of elastic compatibility constraints in matching 
different director and shear modes along the helix. If 
the director at a position z rotates towards the x-axis, it 
is known that the elongation A xx , contraction A yy , and 
shear A xy , are precisely determined by the initial orien- 
tation, (j)Q , and the rotation from it, if the process is to 
be soft. The next slab of material, at z + dz, has the 
initial orientation <po + qqcIz and a different set of soft 
strains A must arise. Material points at y translate to 
Xyy(z) • y and A yy (z + dz) ■ y in the two neighbouring 
slabs along the helix, that is they differ by a relative dis- 
placement {d\yy/ dz)dz ■ y. There is thus a generated 
shear A yz = (d\ yy / dz) ■ y that diverges as the linear y- 
dimension of the sample. We accordingly conclude that 
the transverse contractions are uniform. Such deforma- 
tions, e.g. A zz and A yy , are equal in the first approxima- 
tion, in spite of the apparent anisotropy along the pitch 
axis z (each deformation generates the same energy fceT 
per chain independently of its alignment with respect to 
the director) and are thus both equal to 1/ vA. 

(i) Transverse elongation A xx = A. We consider 
the strain tensor in the following form: 

/AO \ 

A= o i/Va . (4) 
V o o i/Va / 

Although one expects the director rotation in the az- 
imuthal plane (cf. Fig. [l]), there are no associated shear 
strains. Such shears, A xy (z) and A yx (z), would both 
lead to elastic compatibility problems (e.g. a gener- 
ated shear A xz oc y, the sample dimension) and we as- 
sume they are suppressed. The shears A xz (z) and X xz (z) 
are not subject to compatibility requirements. How- 
ever, they should not appear on symmetry grounds, 
which is easily confirmed by direct minimisation. Now 
no = {cos 4>q , sin </>o, 0} and the rotated director after 
deformation is n = {cos <f>, sin 0, 0}. Note that the helix 
is </>o = qoz in the initial undistorted material. After de- 
formation, because of the affine contraction A zz = I/VA, 
the material frame shrinks: z — > z/yX. This has an effect 
of affine contraction of the helical pitch so that q = V\qa 
in all expressions below. With the £ Q and £ implied by 
these n and n, the free energy density Eq. (0) yields 

F± = \* (A 2 + ~ + "-^ [A 2 (r cl - c 2 S 2 ) (5) 

+ i(r c 2 s„ - s 2 c„) - 2VX(r - l)s c s c J , 

where cq and sq are shorthand for cos <j>o and sin tfio; 
analogously, c and s stand for cos <p an d sin <p. The ap- 
pearance of terms linear and quadratic in <f> (or rather 
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sin0 because all values of the azimuthal angle will be 
found along the cholesteric helix) indicate that rotations 
can always lower the energy for A ^ 1. Minimization of 
F± with respect to results in the expression for the lo- 
cal director angle <fr(z) at a given imposed extension A, 
depending on the phase of cholesteric helix: 



tan 2(f> - 



2\ 3 / 2 {r- l)sin 2qz 



(r - 1)(A 3 + 1) cos 2qz + (r + 1)(A 3 - 1) 



(6) 



Initially, all directors at < qz < tt/2 are induced to 
rotate "backward" towards <f> — 0, and all directors at 
—tt/2 < qz < rotate "forward" towards <fi — tt, as the 
imposed deformation A increases, see Fig. H Although 
<f> = and 7r describe equivalent directors, the twist wall 
between these two states becomes more and more sharp. 
Due to the helix imprinting, the orientations <p = at 
qz = and <j) = tt at qz — tt are pinned, as is the middle- 
point of the twist wall at qz = tt/2. As a result, no 
change of the helical pitch can occur. This is in contrast 
with cholesteric liquid crystals, where in a classical prob- 
lem of helix unwinding by electric or magnetic field one 
finds an increase in cholesteric pitch |||| along with the 
coarsening of the helix. 




0.79 1.57 2.36 3.14 

Helix phase qz (radian) 

FIG. 2. The director angle <f) against the cholesteric helix 
phase qoz for increasing strain A = 1 (open circle), 1.15 (trian- 
gle), 1.23 (open square) , 1.25 (shaded square) and 1.5 (shaded 
circle). The polymer chain anisotropy is r = 1.9 and thus the 
critical strain is A c = r 1 ^ 3 ~ 1.24 (see text). At A > A c the 
director pinning at (f> — tt/2 breaks down and a discontinu- 
ous transition occurs, after which the director continuously 
rotates towards the final uniform <f> = 0. 

Examining the Eq.(^) one finds that the denomina- 
tor changes sign and remains negative in the region of 
the twist wall, centered at qz = tt/2 between the values 
(f> = tt/4 and 37r/4. The width of such a wall is 



-A 3 



w 



q V (r-l)(A 3 + l)' 



As the increasing applied strain reaches a critical value 
A c = r 1 / 3 , the wall width w — > and the discontinu- 
ous transition occurs. The director in the mid-point of 
the wall breaks away from the pinning and jumps from 
<j> = tt/2 to (f> = 0, along the strain axis, thus remov- 
ing the topologically constrained twist wall. From this 
point there is no barrier for director rotation towards 



the final uniform orientation with <p = 0, as the last 
two curves in Fig. indicate. The bifurcation occurs 
between the family of <f>(z) curves at A < A c and those 
at A > A c . They are separated by the critical curve <j) c 

given by tan2</> c = —2 tan qz. 

A discontinuous director jump at a critical strain has 
been predicted, and indeed observed in nematic elas- 
tomers stretched at 90° to their initial director no |7|,[l0) ■ 
It was also found that even small deviations from the 
exact perpendicular geometry, or the possibility of soft 
shears in stripe domains lead to a continuous director 
rotation. However, in the present cholesteric case the 
shears are prohibited by compatibility constraints and 
one always finds an exact phase angle cf> — tt/2 along the 
helix - where the center of narrowing twist wall becomes 
pinned from both sides. It is this point that experiences 
a discontinuous jump in order to allow the system to pro- 
ceed to its final uniform equilibrium state. 

(ii) Stretching along the pitch axis A zz = A. We 
now take (cf. Fig. [I]) 



1/VX A xz 
1/VA A yz 
A 



(7) 



No compatibility problem with shears A xz (z) and A yz (z) 
arises from their variation with z along the helical pitch. 
By contrast, their conjugate strains A zx and A zy , which 
would also have to vary with z, would lead to a serious 
compatibility mismatch, e.g. dX zx /dz = d\ zz /dx. We 
therefore assume A zx and A zy are suppressed even though 
in other settings these are the generators of soft elastic 
response. 
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FIG. 3. (a) The geometry of director rotation in response 
to stretching A zz along the helix axis. (b) The angle 6 of 
director tilt plotted against the imposed strain A, Eq. ([ll]) for 
r — 1.3 (triangle) and r = 1.9 (circle). Strain varies from 1 to 
A2 = r 2//3 at which point the alignment is 9 = tt/2, uniformly 
along the former pitch axis. 

In this geometry one expects the director rotation 
out of the azimuthal xy plane, see Fig. ||(a). The ini- 
tial director is, as before, no = {cos qoz, sin q^z, 0}, 
while after deformation the rotated director is n = 
{cos# cos qz, cos9 sin qz, sin#}. As in the case (i), all 
physical dimensions in the deformed sample are scaled 
by the affine strain. In particular, here z — > Xz, resulting 
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in the corresponding expansion of the cholesteric pitch: 
q = qo/X. With the £ Q and I defined by the axes no 
and n, the free energy density Eq. @ now becomes a 
function of three variables: the director tilt angle 9 and 
the two shear strains X xz (z) and A yz (z) (we continue to 
neglect the effects of director gradients and Frank elas- 
ticity). Algebraic minimisation over these components of 
strain tensor is not difficult and results in 



= A- 



(r - 1) sin 29 
1 - (r- l)cos 26» 



cos qz 
sin qz 



(8) 



in phase with the azimuthal angle along the helical pitch. 
Eq. (|^) describes small distortions in the xy plane, per- 
pendicular to the helix axis, rotating following the initial 
orientation no. On substitution of these optimal shears 
back into the free energy density one obtains 



Fu 



1 



A 2 



2 + (r - 1) sin 2 



2 r " Vl + (r- 1) sin 2 9 
Fn expands at small tilt angle 9 as 
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F„ « ^(A 2 + 2/A) - \^{r 1)(A 2 - 1/A), 



(9) 



(10) 



that is the director starts to rotate down to define a cone 
of semiangle tt/2 — 9 immediately as the strain A > 1 
is imposed. The equilibrium director tilt is obtained by 
minimisation of the full free energy density Fu (9): 



sin a = 



A 3 / 2 - 1 _ 
- 1 



• /A 3 / 2 -l , s 
arcsm-w — . (11) 



The director rotation starts and ends in a characteristi- 
cally singular fashion Fig. [|(b) (reminiscent of the univer- 
sal opto-mechanical response seen in nematic elastomers 
PH). The rotation is complete with the director aligned 
along the extension axis (9 = tt/2) at A = r 2 / 3 which, for 
some elastomers, can be a very large extension. 

In contrast to conventional cholesteric liquid crystals, 
we have altogether ignored effects of Frank elastic energy. 
The most compelling evidence for this is the very stabil- 
ity of the imprinted helical state in the face of the Frank 
penalty |i^2?o- The argument for this relies upon the 
great difference in characteristic length scales, the elas- 
tomer penetration depth, more accurately expressed as 
£ ~ -^jy/K/fi [cf. Eq. (||)], and the director modulation 
wavelength estimated by the helical pitch p = 7r/g 3* £• 
There are two possibilities to alter this inequality - by 
increasing the penetration depth £ (either by making a 
weaker gel, or a less anisotropic one), or by locally in- 
creasing the director gradient (for instance, in the ever 
narrowing twist wall, Fig. ||) . 

One can estimate how weak a gel must be for the Frank 
elasticity to intervene in our analysis. When £ ~ p, for 
example with a pitch p ~ 4 x 10 -7 m, then a rubber modu- 
lus of only \i ~ 60 J/m 3 is required (assuming [r— 1] ~ 1). 
Nematic elastomers typically have /i ~ 10 3 — 10 5 J/m 



and their cholesteric analogues would clearly find Frank- 
elastic effects minor. However, an elastomer with a rea- 
sonable /i ~ 10 3 J/m 3 would feel the director gradients 
when its polymer chain anisotropy becomes as low as 
r = ~ 1.25. Such a value is easily reached in side- 

chain liquid crystal polymers, especially near the clearing 
point §. 

Another interesting test of the role of Frank elastic- 
ity is in the twist wall described in our case (i). The 
width of the wall decreases to zero, and the Frank en- 
ergy density grows, being maximal at the centre of the 
wall. There it is ^K 2 qoX 3 (r — l) 2 /(r — A 3 ) and diverges 
at the critical strain A c . Therefore, the local analysis 
of Eqs. (||)-([|) is on ly valid outside the region of strain 
AA pa ^ 3 r 2/^ j (loO 2 around A c . In a typical hard ne- 
matic rubber this is a narrow region of AA < 0.01, but in 
a weak gel with low chain anisotropy it may become more 
substantial. Moreover, the finite width of the twist wall, 
demanded by the Frank gradient energy, raises the ques- 
tion of topological mechanism for eliminating the twist 
stored in the cholesteric helix, perhaps by a disclination 
loop expansion. 

To summarise, we have predicted a qualitatively new 
response of an elastomer with chiral cholesteric mi- 
crostructure to applied fields that is different from clas- 
sical cholesteric liquids. Likewise, the chiral imprinting 
and its modification by elastic fields is a new effect in rub- 
bers and solids. One could envisage tuning these effects 
by the use of solvents (with or without chiral power) and 
by other fields affecting the director, for instance electric. 
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